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ABSTRACT 

We develop a simple star formation model whose goal is to interpret the emerging body of observational 
data on star-forming galaxies at z > 6. The efficiency and duty cycle of the star formation activity within 
dark matter halos are determined by fitting the luminosity functions of Lya emitter and Lyman-break galaxies 
at redshifts z — 5 - 6. Our error budget takes proper account of the uncertainty arising from both the spatial 
clustering of galaxies and the commonly-used Poisson contribution. For a given survey geometry, we find a 
cross-over luminosity below which clustering provides the dominant contribution to this variance. Using our 
model parameters we predict the likely abundance of star forming galaxies at earlier epochs and compare these 
to the emerging data in the redshift interval 7 < z < 10. We find that the abundance of luminous Lyman-break 
galaxies in the 500 Myr between z — 6 and 10 can be naturally explained by the hierarchical assembly of dark 
matter haloes; there is only marginal evidence for strong physical evolution. In contrast, the first estimates of 
the abundance of less luminous star forming galaxies at z = 9- 10 are higher than predicted and, if verified by 
further data, may suggest a top-heavy stellar mass function at these early epochs. Although these abundances 
remain uncertain because of the difficulty of spectroscopic confirmation and cosmic variance, even a modest 
improvement in survey capability with present or upcoming facilities should yield great progress. In this 
context, we use our model to consider those observational techniques that hold the most promise and make 
predictions for specific surveys that are, or will soon be, underway. We conclude that narrowband Lya emitter 
surveys should be efficient on searches at z — 7-8; however, such conventional surveys are unlikely to detect 
sufficient galaxies at z — 10 to provide useful constraints. For this reason, gravitational lensing offers the best 
prospect for probing the z — 10 universe prior to JWST 

Subject headings: cosmology: theory — galaxies:formation — galaxies: high-redshift 



1. INTRODUCTION 

The recent discovery of star-forming galaxies at high red- 
shifts, z > 6, represents an emerging frontier in our under- 
standing of the early stages of galaxy formation. Such stud- 
ies aim to determine the role that young galaxies may play 
in completing cosmic reionization, as well as to define more 
clearly how feedback and other poorly -understood processes 
shape the early distribution of galaxy luminositi es and sizes 
from whic h later systems evolve (see reviews by|LoebJ[2006; 
lEUTsll2007h . 

Considerable observational progress is being made through 
ambitious campaigns undertaken with the Hubble and Spitzer 
Space Telescopes and large ground-based telescopes. As 
a result, luminosity functions and stellar mass distribu- 
tions are now available for various star-forming popula- 
tions observed at z — 5 - 6. These include the continuum 
'drop - outs' or Lyman break galaxies (L BGs.Bouwen s etaLl 
l200a lYan et all 120061 lEvles et all [20061: IStark et al l 12006), 
and the Lyman-alpha em i tters (LAEs. | Santos et al.l l2004t 
Malhotra & Rhoadsl l2004t iHu et al.l bOOiT Tshimasa ku et al l 
2006; Kas hikawa et al.l 120061) . located either spectroscopi- 
cally or via narrow band imaging. Alongside these achieve- 
ments, the first constraints are now emerging on the abun- 
dance of equivalent systems at 7 < z < 10 (Richarcfet al.1 



20061; IWillis & Courbinll2005l: iBpuwens & Il lingworthll2006l: 
Ive et al.ll2006l;ICubv et a l. 2006; Stark et al. 200J) 



Several questions arise as observers continue to make 



MS 



' Department of Astrophysics, California Institute of Technology, 
105-24, Pasadena, CA 91125; dps@astro.caltech.edu 

2 Astronomy Department, Harvard University, 60 Garden Street, Cam- 
bridge, Massachusetts 02138, USA 



progress. First, what is the physical relationship between 
the Lyman Break Galaxy (LBG) and Lyman-alpha emitter 
(LAE) populations? This is important in interpreting the 
quite significant differences that have been observed between 
the properties of the two classes. Second, what redshift 
trends are expected for these populations? Some authors 
(e.g. iKashikawa eTaTI 120061: feouwens & Illingworthl 120061) 
have claimed strong evolution in the abundances with redshift. 
However, in the absence of any theoretical framework it is dif- 
ficult to assess the significance of these claims. Finally, given 
the abundance of star forming galaxies at z — 5 - 6, what is 
expected at z — 7- 10, the current observational target? And 
what is the optimum observational strategy for finding those 
sources which will be valuable in constraining the epoch of 
cosmic reionization and the properties of young galaxies? 

The present paper is motivated by the need to answer 
these questions given the improving observational situation 
at z ~5-6. As the time interval between z=6 and z=10 is short 
by cosmic standard s (~470 Myr in the WMAP3 cosmology, 
ISpergel et al.l2006l) . the growth of the halo mass function over 
this redshift range is well understood. Accordingly, it is prac- 
tical to fit the z — 5 - 6 observations in the context of a simple 
star-formation model, thereby deducing the likely differences 
between the LBG and LAE populations, and to then use such 
an empirically calibrated model to predict their likely abun- 
dances at z — 10. 

The model we adopt assumes all early star-forming galax- 
ies observed during this relatively short period are primarily 
triggered into action by halo mergers. The key parameters 
governing their visibility are thus the efficiency of star forma- 
tion and the duty cycle of its activity. An observed luminosity 
function at z — 5-6 thus provides a joint constraint on the 
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star formation efficiency (which determines the rate at which 
baryons are converted into stars) and the duty cycle of activity 
(which determines the fraction of halos occupied by visible 
galaxies). We consider such a simple model to be comple- 
mentary and perhaps more intuitive than full ab initio numeri- 
cal simulations which must also assume sub-grid presriptions 
for st ar formation (e.glNagamine et al.ll2005t iGnedin & Fan! 
[2005 iFinlator et alJl2006h . The goal is to infer the likely red- 
shift trends in the context of emerging data and to use the 
model to evaluate the future observational prospects, partic- 
ularly in the optimal design of surveys to locate sources at 
Z — 10 or so. 

Several authors have alread y made good p r ogress with 
such semi-analytic models. Le Delliou et al. (20051) and 
Diikst raet alj ([2006 ) have attempted to fit the luminosity dis- 
tribution of LAEs. Dijkstra et al. argue that the evolution ob- 
served by Kashikaw a et al.l J2006) between z =5.7 and z=6.5, 
claimed as arising from changes in the intergalactic medium, 
may be underst ood instead through simple growth in the halo 
mass function. Ma o et all (120061) have extended the method 
to contrast the properties of LAEs and LBGs; they find LBGs 
reside in a wide range of halo masses (10 10 to 10 12 M Q ), 
whereas LAEs reside within a narrower range (~ 1O U M0). A 
key inference from their model is t he short duty cycle of ac- 
tivity in the most luminous sources. ISamui et all (120061) have 
argued that the evolution observed at 6 < z < 10 in the Hubble 
Ultra Deep Field (UDF) can be attributed to evolution in the 
underlying dark matter halo number density without requiring 
a dramatic change in the nature of star formation; in contrast, 
they find that the large abundance o f z — 9 LBGs discove red 
in the gravitational lensing survey of Richard ^t al.l (120061) re- 
quires significant evolution in the stellar initial mass function, 
the reddening correction, and the mode of the star formation. 

The present paper continues the earlier work. We focus not 
only on explaining the growing body of data at z — 5 - 6 but 
will also include the emerging data at higher redshift to see 
how it agrees with our model predictions. A crucial ingredient 
in finding the model parameters that best fit the observational 
data is the error budget on that data. The commonly-used scat- 
ter due to Poisson fluctuations must be supplemented by the 
spatial clustering of galaxy halos within the volume of each 
survey. Here we will provide new error bars for existing data 
and show that the cosmic variance due to clustering dominates 
over Poisson fluctuations below a particular galaxy luminos- 
ity for a given survey geometry. This understanding can be 
used to optimize the flux sensitivity and strategies of future 
high-z surveys. 

The plan of the paper is as follows. In §2 we introduce 
the basic ingredients of our model for star forming galaxies 
at high-redshifts. In §3 we calculate the error budget for ob- 
servations of LBGs and LAEs at high redshift. The model is 
calibrated against existing data for LBGs and LAEs at z ~5-6 
in §4 and used to discuss the emerging data at earlier epochs in 
§5. In §6 we discuss the implications of our model for further 
observational campaigns with current and projected facilities. 
In §7, we summarize our conclusions. 

Throughout the paper, we have assumed a flat universe 
andjT2 m ,f2A)=(0.24,0.76) following results in Soer gel etail 
(2006t). All magnitudes are given in the AB system. 

2. A PHYSICAL MODEL FOR HIGH REDSHIFT STAR FORMING 
GALAXIES 

The rationale of this paper is to empirically-calibrate, using 
the data now available at z — 5 -6, the parameters of a simple 



model that describes the evolving luminosity function of star- 
forming galaxies (both LBGs and LAEs) over the quite short 
time interval corresponding to the redshift range 5< z <10. 
Such a model can then be used to make predictions for the 
upcoming z — 7- 10 surveys. 

The model assumes star formation at these early epochs 
is primarily triggered by the well-understood rate at which 
dark matter halos coalesce. We assume that the ratio of 
baryons to total mass in halos above some minimum mass 
dWvithe & Loebl 120061) follows the universal value flt,/Q m . 
Baryons are subseque ntly converted to sta rs w ith an efficiency 
given b y f*. Following lLoeb et ail d2005) and lWyithe & Loebl 
(2006), we define the star formation timescale, Ilt, as the 
product of the star formation duty cycle, eoc, and the cosmic 
time, te = 2/3H at that redshift. Using these ingredients, the 
star formation rate M* is related to halo mass Mh a i as follows 

xa \ f * X ( fi b/^ m ) X M halo 

M + (M ha i ) = '- . ( 1 ) 

tLT 

For comparison to LBG samples, the star formation rate de- 
fined above is converted to the luminosity per unit frequency 
at 150 A following the prescription presented in Mad au et alj 
(fT99l : L 150 o = 8.0 x 10 27 (M*/M Q yr 1 ) erg s^Hz" 1 . This 
conversion factor assumes a Salpeter initial mass function 
(IMF) of stars; if the IMF is more top-heavy than the Salpeter 
IMF, the far-ultraviolet luminosity will be greater for a given 
M*. 

Comparison to LAE samples requires converting the star 
formation rates derived above to a Lya luminosity. We do 
this assuming that two-thirds of all recombining photons re- 
sult in the emission of a Lya photon (case B recombination). 
The ionizing photon production rate is calculated from the 
star formation rate for a given IMF and metallicity. We fix 
the metallicity at 1/20 solar and assume a Salpeter IMF which 
yields N 7 = 4 x 10 53 i onizing photon s per second per star for- 
mation rate in M Q / yr ( Schaerer 2003). An extreme top-heavy 
Population III IMF would produce 3 x 10 54 ionizing photons 
dBromm et al.l 120011; ISchaereri 120031) . and we will consider 
such an IMF in later sections. Since the Lya photons are as- 
sumed to be produced via recombinations, only the fraction 
of ionizing photons which do not escape into the intergalactic 
medium, (1 -f; p ), produce Lya photons. Furthermore, only a 
fraction, T^ ya , of the emitted Lya photons escape the galaxy 
and are transmitted through the intergalactic medium (IGM). 
With this prescription the Lya luminosity is related to the halo 
mass as follows: 

2 

Lhya = ^hv L y a N 7 T hya (l -/ ip )M*. (2) 

A substantial change in the IGM transmission parameter T Lya 
is expected to signal the end of reionization. 

We also consider a more advanced model, incorporating 
the effects of supernova feedback on the luminosity func- 
tion of star-forming galaxies. Since supernova feedback can 
significantly reduce the efficiency of star formation in low- 
luminosity galaxies, it is particularly important to consider 
when predicting the efficiency of future surveys for high- 
redshift galaxies aimed at detecting intrinsically fainter sys- 
tems, rather than the luminous and rare objects that have been 
detected thus far, Fol lowing the scaling relations presented in 
lDekel& Wool ([2003). we assume there is a critical halo mass 
at each redshift below which the star formation efficiency be- 
gins to decrease due to feedback. The star formation effi- 
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ciency (which we now call rj(M)) is a function of halo mass, 
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where Mhaio,crit represents the critical value. At a given red- 
shift, the critical halo mass can be related to a critical halo 
velocity. In the local universe, observations sugges t a critical 
halo velocity of - 100 km s" 1 (iDekel & Wool2003l) : we adopt 
this for our high-redshift models, assuming that the physics of 
supernova feedback depends only on the depth of the gravita- 
tional potential well of the halos. 

Finally, in §6 we also consider the expected evolution of 
galaxy sizes. This is particularly important for gauging the ef- 
ficiency of future surveys utilizing adaptive optics in detecting 
galaxies at z > 7-20. If a galaxy is resolved in a particular ob- 
servation, the signal-to-noise ratio for the detection increases 
as the size of galaxy decreases since a smaller aperture (with 
less noise) is required to encircle its flux. Next-generation 
adaptive optics systems on thirty-meter class ground-based 
telescopes will offer a resolution of >9 milliarcseconds at 
1 . 1 /im, corresponding to ~50 pc at z —7. With such exquisite 
angular resolution, the size of early galaxies is likely to be a 
limiting factor in their detection. 

We assume that the extent of the stellar region of a galaxy 
at a particular epoch is a constant fraction of the size of the 
host dark matter halo. The virial radius of a halo is given by 



GM,,„. 
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(4) 



where M V! > is the halo virial mass and H(z) = Ho[f2A,o + (l - 
f2 A . -^o)(l +z) 2 +^o(l +z) 3 ] l l 2 is the Hubble canstant at red- 
shift z- At very high redshift, the halo virial radius scales as 
(1+z)" 1 for fixed halo mass. Follow ing the parameterization 
presented in Barkana "& Loebl(l2000h . a typical galaxy brighter 
than 1 nJy between 5 < z < 10 will have a disk radius of 0"1- 
Of'2, with the range depending on the efficiency with which 
baryons are converted to stars. 

Observations lend support to this simple scaling of disk ra- 
dius with redshift: z — 2 - 6 dropouts in the UDF and Hubble 
Ultra Deep Field Parallels (UDF- P) are best fit by a (1 +zT 1 
power law for a fixed luminosity iBouwens et alj d2004l) . The 
mean half-light radius of 0.3-1.0 L* >z =3 LBGs at z — 6 is 0.8 
kpc, corresp onding to P."! at z ^ 6, co nsistent with the model 
presented in Barkan a~& Loebl d2000l) . To make predictions 
for future observations in §6, we thus assume the mean size 
of galaxies of similar luminosity scales as Of 1(1 +z/7) _1 . 

3. THE EFFECT OF VARIANCE IN DEEP SURVEYS 

Interpreting and planning observations of galaxies in the 
reionization era requires an accurate understanding of uncer- 
tainties arising from both Poisson errors and fluctuations in 
the large-scale distribution of galaxies. For example, one 
of the main motivations for galaxy surveys at z ~7-10 is 
the question of whether the impr int of reionization can be 
seen in the evolutio n of LAEs (iMalhotra & Rhoads l |200l 
Dijkstra et al.1 l2006h or d warf galaxies dBarkana & Loel 
2006tlBabich & Loebll2006l) . However, to answer these ques- 
tions, it is important to ensure that the variance is less than the 
claimed evolution. In this section, we develop the formalism 
necessary to compute the variance for narrowband, dropout, 
and spectroscopic lensing surveys, and we then apply this for- 
malism to recent surveys. An analysis of clustering variance 



was presented in ISomerville et all (120041) . We improve upon 
two simplifications made in that work, both of which we dis- 
cuss below. 

The variance due to Poisson errors is given by A/), where 
Mi is the number of galaxies in luminosity bin i. To com- 
pute this we determine the number density of galaxies (either 
LBGs or LAEs) as a function of luminosity using the best- 
fitting models described in §4, from which we can determine 
the predicted counts for a given survey volume. When charac- 
terizing the variance at redshifts where data is sparse (z > 7), 
we assume no evolution in the model parameters from z — 6. 
This method clearly has its limitations when probing new pa- 
rameter space (e.g. low luminosities or high redshifts) where 
the abundance of galaxies is not well-known. However, given 
the lack of data, we consider it to be the simplest approach. 

The clustering of galaxies in overdense regions causes fluc- 
tuations in galaxy counts, often referred to as cosmic vari- 
ance. Determining this variance requires knowledge of the 
mass of the dark matter haloes that host the observed galaxies. 
If the correlation function of the galaxy population is known 
the clustering variance can be predicted for a given field of 
view. Moreover, if the data set is sufficiently large, the corre- 
lation function can be derived separately for bright and faint 
galaxies, thereby showing how the dark matter halo mass (and 
hence the clustering variance) varies with galaxy luminosity. 
A detailed spatial correlation function analysis is very chal- 
lenging for current observations at z > 6; hence an alternative 
method is needed for deriving the clustering variance. 

Using the LBG and LAE model described in §2, the halo 
mass can be determined as a function of galaxy luminosity. 
The clustering variance can then be calculated for a given halo 
mass by taking the product of the variance of dark matter in 
the survey volume and the bias factor associated with halos of 
a given mass M. With this method, we compute the clustering 
variance as a function of galaxy luminosity for narrowband, 
dropout, and spectroscopic lensing surveys, taking into con- 
sideration the survey geometry specific to each survey. 

First, we compute the variance of dark matter in a given 
smoothing window as follows: 



cr 2 (r)= / P(k)W 2 (k)d 3 k 



(5) 



where P(k) is the power spectrum of density fluctuations ex- 
trapolated to z=0 as a function of wavenumber k and W 2 (k) is 
the Fourier transform of the window function in real space. 
The form of window function depends on the survey ge- 
ometry; in the following subsections, we detail the window 
functions adopted in our analysis. Non-linear corrections 
to the power spectrum and probability distribution become 
important if the variance is larger than unity. In comput- 
ing the variance, we adopt the non-linear power spectra us- 
ing the halo model fitting functions presented in Smi th et alJ 
(2003). While fluctuations are in the linear regime, their 
probability distribution is Gaussian. However, in the non- 
linear regime, the probability distribution functi on appears 
to be well-described by a log-normal distribution ( Kay o et ail 
1200 lb . Confidence intervals are thus determined via the geo- 
metric mean and standard deviation; the one-sigma upper and 
lower bound are given by exp(/.t + a) and exp(/.t - a), respec- 
tively, where /i and a are the mean and standard deviation of 
the logarithm of a given density Auction, In S. If the standard 
deviation of the log-normal distribution is much less than one, 
&ln *C 1, the probaility distribution function reduces to Gaus- 
sian with confidence intervals given by // ± a. 
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FIG. 1. — Variance in a narrowband Lya survey. The total variance (solid 
line) is the sum of the variance from Poisson noise (dotted line) and clustering 
fluctuations (dashed line). The survey specifications adopte d in the left panel 
are equivalent to Subaru survey for Lya emitters at z = 5.7 (Shima saku et all 
2006). The clustering fluctuations in such a narrowband survey are very small 
(less than 6%) in the luminosity range over which Lya emitters are detectable. 



The clustering variance in the distribution of galaxies is 
then estimated by multiplying the variance in dark matter by 
the halo bias, which is defined as the ratio of the rms fluc- 
tuations of haloes to that of dark matter. We adopt the halo 
bias formu l a deri ved for the ellipsoidal collapse model by 
IShefhetal.l(l2Q0l . 

A key assumption in the clustering variance formalism de- 
scribed above is that there is not more than one galaxy per 
halo. More complex occupation numbers are possible, but 
since cooling is efficient at such high redshifts, the simplest 
case is that with one galaxy per halo. 

Now we introdu ce two key differences between our ap- 
proach and that of Som erville et al.1 (120041). First, the sur- 
vey ge ometry was assumed to be spherical in lSomerville et all 
(2004). However, different survey geometries may have sub - 
stantially different power spectra (Kai ser & Peacocklll991l) . 
This is a particular concern for strong lensing s urveys that 
utilize longslit spectroscopy ( Santos et al.l l2004; Sta rk et al.l 
120071) . Second, the observed number density of the popu- 
lation was assumed to be equivalen t to the number density 
of the underlying dark matter halos (ISomerville et al.ir2004l) . 
This assumption could be in error if the star formation duty 
cycle, £dc, is significantly smaller than unity. In this case, the 
observed number density of galaxies would be less than that 
of their host dark matter haloes by a factor of eoc- Depending 
on the slope of the mass function, this would overestimate or 
underestimate the cosmic variance. 

In the following subsections, we apply the formalism de- 
scribed above to recent surveys for LBGs and LAEs. For each 
survey, we compute the typical amplitude of clustering and 
Poisson fluctuations in the relevant observational regime. In 
addition, we determine the cross-over luminosity, L e , defined 
as the luminosity above which the Poisson error dominates 
that from clustering fluctuations. 

3.1. Variance in Narrowband Surveys 

The geometry of a narrowband LAE can be approximated 
as a rectangular parallelepiped with a comoving volume of 
a x a y a z , where a x and a y are the comoving distances corre- 
sponding to the areal field of view of the survey and a z is the 
comoving line-of-sight distance of the survey. The window 
function in k-space is simply the Fourier transform of a rect- 
angular top-hat with dimensions corresponding to the survey 



sin(k x a x /2) sm(k y a y /2) sin(k z a z /2) 
(k x a x /2) (My/2) (k z a z /2) 



(6) 



Using the window function defined above, we evaluate the 
expected variance in narrowband surveys for LAEs at z — 6. 
We focus our analysis on the narrowband surveys conducted 
at z = 5.7 and z = 6.5 in the Sub aru Deep Field (SDF) . Using 
the wide-format, Suprime-Cam (Mivaz aki et alj [2002) on the 
Subaru Telescope, the observations sample an area of 34' x 27' 
in the NB816 and NB921 narrowband filters. The central 
wavelength and FWHM of the NB816 and NB921 filters are 
(8150 A, 120 A) and (9196 A, 132 A) respectively. 

In Figure Q] we plot the variance in the z = 5.7 narrowband 
survey of the SDF. The Lya luminosities probed in the Subaru 
survey range from ~ 10 — 10 43 erg s" 1 . The clustering vari- 
ance is less than ~0.01 over this luminosity range, resulting 
in less than 10% uncertainty in the observed counts. Poisson 
errors dominate the clustering errors for sources brighter than 
the cross-over luminosity of 10 41 ' 7 erg s" 1 (Table 1). Since this 
is intrinsically fainter than the luminosity limit of the SDF, 
this survey is dominated by Poisson errors. 

3.2. Variance in Lyman-break Surveys 

We focus our dropout survey analysis on the two 
16'xlO' Great Observatories Origins Deep Surveys 
(GOODS) of the Hubble Deep Field North (HDF-N) and 
Chandra Deep Field South (CDF-S) and the 3.4' x 3.4' Hub- 
ble Ultra Deep Field (UDF). The redshift distribution of 
dropouts depends on the filters and color-cuts used in their se- 
lection. While the typical color-selection criteri a for i— drops 
select galaxies between z = 5.5 and z = 7.0 ( Bun ker et al.l 
l2004al) . the effective distance sampled along the line-of-sight 
is less than the total comoving radial distance in this redshift 
interval because of incompleteness arising from objects be- 
ing scattered faintward of the magnitude limit or out of the 
color-selection window. This incomplet eness has been quan- 
tified in both GOODS and t he UDF dBunker et all l2004at 
iBouwens & Illin gworth 2006), allowing an effective volume 
to be derived for each survey. We approximate the geome- 
try of the LBG survey as a rectangular parallelepiped with a 
characteristic line-of-sight distance equal to the ratio of the ef- 
fective volume and the survey area and with dimensions in the 
plane of sky corresponding to the field-of-view of the survey. 

The variance of the GOODS and UDF surveys for z — 6 
LBGs is presented in Figure [2] The GOODS survey is sensi- 
tive to sources brighter than 5 x 10 28 erg s" 1 Hz" 1 (correspond- 
ing to a lOcr limit of z^ B =26.5). At this limit, the clustering 
variance is only 0.003, and the total variance is dominated 
by the Poisson term. The clus tering fluctuations are s lightly 
lower than those estimated in ISomerville et al.l (l2004h . due 
largely to the more realistic source geometry. In the UDF, 
both the clustering and Poisson error are larger due to the 
smaller survey volume. The 8er li miting far UV lum inosity 
in the UDF is 8 x 10 27 erg s" 1 Hz" 1 dBunker et al.ll2004bl) . The 
total variance near this limit is dominated by clustering fluc- 
tuations, which contribute an uncertainty of ^15-20% to the 
observed densities. 

The field-t o-field fluctuations of z — 6 LBGs have been 
measured in IBouwens et al.l d2006h by degrading the depth 
of the UDF to that of the two UDF-parallel fields and subse- 
quently comparing the number of i-dropouts in each field. The 
density of i-dropouts selected in the (degraded) UDF is similar 
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FIG. 2. — Variance in dropout LBG survey. The total variance (solid 
line) is the sum of the variance from Poisson noise (dotted line) and clus- 
tering fluctuations (dashed line). In the case of a survey such as the HST- 
GOODS observations of the CDF-S and HDF-N (limiting luminosity of 
5 X 10 2s ergs~'Hz -1 ), the clustering fluctuations of z — 6 LBGs are greater 
than ~6% (top panel). In the Hubble Ultra Deep field (limiting luminos- 
ity of 8 X 10 27 ergs~'Hz~'), clustering fluctuations result in a slightly higher 
uncertainty of >15-20% (bottom panel). 

to that in the first parallel (50.2 and 42.6, respectively) but is 
significantly greater than that in the second parallel (27.8 vs. 
11.4). The latter comparison implies field-to-field variations 
on 7 arcmin 2 at the magnitude limit (zao = 28.6 at 8cj) of the 
second UDF parallel are ~ 40%. The observed field-to-field 
variations result from both Poisson (19% for the degraded 
UDF and 30% for the second UDF parallel) and clustering 
fluctuations (19-26% over 7 arcmin 2 ). When these uncertain- 
ties are accounted for, the two measurements are consistent at 
the 1.3cr level. 

3.3. Variance in Lensed Longslit Spectroscopic Surveys 

The geometry of a longslit spectroscopic survey can be ap- 
proximated as a rectangular parallelepiped. However, in the 
case of a strong lensing longslit spectroscopic survey, the 
geometry is potentially slightly more complex. For a lens- 
ing survey, the slit geometry only corresponds to the image 
plane; however, we are interested in the geometry of the sur- 
vey in the source plane, which can be calculated accurately 
via ray tracing if a reliable mass model is available. While 
the source plane geometry depends on the location of the 
lo ngslit relative to th e critical line, for typical clusters studied 
in lStark et al.l (1200 7) it is well-approximated by a rectangular- 
parallelepiped (J. Richard 2006, private communication). The 
source plane area is reduced by a factor of the lensing magni- 
fication, M.\ further, the magnification is not isotropic and is 
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FIG. 3 . — Variance in a z = 10 lensed longslit spectroscopic survey. The to- 
tal variance (solid line) is the sum of the variance from Poisson noise (dotted 
line) and clustering fluctuations (dashed line). The uncertainty from cluster- 
ing fluctuations is significantly greater than the Poisson noise for all lumi- 
nosities and is nearly 100% for sources with Lya luminosities of 10 40 -10 41 
and erg s _1 . 

strongest perpendicular to the cluster critical line. Hence, as- 
suming the longslit is oriented along the cluster critical line, 
the slit-width is compressed more than the slit-length. For 
the computations that follow, we assume the source plane 
slit width, a x , is related to the image plane slit width, a\ 
by a x =a' x / V2JA, and likewise, the source plane slit width is 
given by a v =a' y \/2/ 'M, in a greement with typical slit posi- 
tions from lStark et all (120071) . 

As with the narrowband survey, the appropriate window 
function is a three-dimentional rectangular top-hat in real 
space, which in k-space corresponds to the product of three 
sine functions (see Equation|6]l. 

A typical near- infrared spe ctrometer has dimensions of 
0.76" x 42" dStarketalJ l2007), which at z ~ 9 corresponds 
to a comoving distance of 6 kpcx610 kpc, assuming a me- 
dian magnification factor of Ai=2Q. In each cluster, an area is 
mapped out around the critical line; assuming six slit positions 
are observed, this results in a total survey area of 0.02 Mpc 2 
per cluster. If observations are conducted in the J-band be- 
tweenz= 8.5 andz= 10.4, the comoving line-of-sight distance 
spanned is 479 Mpc. We compute the clustering variance ex- 
pected over the total survey volume for a fifteen cluster survey. 
In the near future, significantly m ore clusters will b ecome 
available for strong lensing surveys dEbeling et al.ll2003l) . 

The clustering fluctuations are significantly larger in the 
spectroscopic lensing survey than in either the traditional nar- 
rowband or dropout surveys due to the much smaller survey 
volumes. While such surveys may offer the only prospect of 
detecting galaxies at z — 10, clearly there will still be large 
uncertainties in their abundance due to cosmic variance. The 
development of larger and more sensitive near-infrared spec- 
trometers is necessary to increase the survey volume obtain- 
able in a reasonable time allocation. 

4. MODEL CALIBRATION USING Z ~ 5 - 6 OBSERVATIONS 

We are now in a position to use our model and our improved 
understanding of the effects of cosmic variance to constrain its 
parameters using observations of the luminosity functions at 
z ~ 5-6. By doing this independently for LBGs and LAEs 
we will hopefully gain valuable insight into the physical dif- 
ferences between these two star-forming populations. 

4.1. Lyman-Break Galaxies 

Lyman-break galaxies are perhaps more straightforward to 
model than Lya emitters because of the complex resonant in- 
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teraction of Lya photons with neutral hydrogen which occurs 
in the latter population. In our model of LBGs described in §2, 
there are two free-parameters: the star formation efficiency /* 
and the duty cycle eoc- First we determine these parameters 
by reproducing the observations at z — 6. We later use this 
model to consider whether the emerging data at z — 7 - 10 
requires any adjustment. Significant evolution in the model 
parameters between the two redshifts might signify some ex- 
ternal phenomenon, such as the reionization of the intergalac- 
tic neutral hydrogen. Alternatively it could cast doubt on the 
reliability of the observations. 

We compute a grid of LBG luminosity functions by vary- 
ing f+ and eoc- The comoving number density of galax- 
ies predicted by the models, n mo( j, is compared to the ob- 
served value, n b s , in each of the N luminosity bins, and a 
likelihood of a given set of parameters is defined such that 

£(f*,e DC ) = exp[-0.5x 2 ], where x 2 = £jli(nobs,i- n mod,i)- The 
1-sigma uncertainty in the observed densities include the con- 
tribution from cosmic variance (see §3) in addition to that 
from Poisson noise. 

We first apply our model to the observed abund ance of 
LBGs at z — 6, as compiled bv lBouwens et al.l d2006l) . In Fig- 
ure |U we show the likelihood contours at 64% and 26% of 
the peak likelihood (corresponding to 1- and 2-er for a Gaus- 
sian distribution). The maximum likelihood and 1-a con- 
fidence intervals are (/*,e DC ) = [0.13^,0.20!^], in rea- 
sonable agreement w ith a similar fit to these observations in 
Wvit he& Loebl(l2006T) . When supernova feedback is allowed 
to decrease the star formation efficiency in low-mass halos 
(see §2 for details), the best-fit parameters change slightly: 
(A,£dc) = [0.16^^,0.25^;,^] and the likelihood increases 
by almost a factor of two. The strong degeneracy between 
the duty cycle and the star formation efficiency arises because 
an increase in the star formation efficiency requires a longer 
star formation timescale (and hence larger star formation duty 
cycle) to produce the same far-UV luminosity for a given halo 
mass. 

Although there is some degeneracy between the best-fitting 
star formation efficiency and duty cycle at z — 6, the range of 
values can be physically understood. A duty cycle of 20% at 
z ~ 6 corresponds to a star formation lifetime of ~ 200 Myr 
which is only slightly larger than the dynamical time of viri- 
alized halos (or the duration of equal-mass mergers) at that 
redshift. Our simple model thus suggests that at z — 6, star for- 
mation is proceeding on roughly the same timescale it takes 
virialized baryons to settle to the center of the galaxy. A star 
formation efficiency of ~ 13% is reassuringly similar to the 
ratio between the global mass density in stars and baryons in 
the present-day Universe (Fukugi ta et al.|[l998"l) . 

An independent check on the inferred duty cycle could con- 
ceivably be obtained if the spectral energy distribution was 
known for a large sample of z — 5 - 6 L BGs. Fitting these 
with a grid of population synthesis models lBruzual & Chariot! 
(2003) allows, in principle, the estimation of the stellar mass, 
dust extinction, and luminosity-weighted age of representa- 
tive galaxies. Unfortunately, the ages inferred via this tech- 
nique have large systematic uncertainties due to the inabil- 
ity of the population synthesis models to constrain the pas t 
star formation history dEvles et al.ll2QQ6t IShaplev et al.ll2005l) . 
Taking the star formation hi stories that minim ize the \ 2 fit to 
the SEDs of LBGs at z ~ 6, lEvles et al.l d2006l) find a median 
age of 500 Myr for those objects detected in the rest-frame 
optical with Spitzer (and hence the most massive objects). A 
stacking analysis of the least massive LBGs in their survey 



1.00 




n.m I , I 

0.01 0.10 1.00 

FIG. 4. — Top: Confidence intervals on star formation efficiency /* and 
duty cycle eoc in a simple theoretical model for the observed abundances of 
LBGs at z=6. The likelihood contours are 64% (blue dashed line) and 26% 
(red solid line) of the peak likelihood. Bottom: Same as above except with 
supernova feedback included in model. 

indicates that these objects have ages of ~60 Myr. Recall- 
ing that the duty cycle is equal to the star formation timescale 
divided by the cosmic time, these inferences suggest that the 
duty cycle lies in the range 6-50% and perhaps increases with 
the mass of the galaxy. 

A further check is provided by limited data on the cluster- 
ing of LBGs. The halo masses probed in the iBouwens et alj 
(2006) compilation in the GOODS and UDF surveys are 
7 x 10 9 Mq-3 x 10 11 M Q according to the simple model 
we have adopted. These values are consistent with cluster- 
ing analysis of z — 6 LBGs in GOODS, whi ch suggest that 
the ho sting dark matter haloes are ~ 10 11 Mq (lOverzier et alJ 
2006). Reionization would be accompanied by a dramatic in- 
crease in the cosmological Jeans mass and and th erefore in 
the minimum galaxy mass (Wvith e & Loebl 120061) . A direct 
detection of this effect requires finding galaxies in dark mat- 
ter haloes over an order of magnitude less massiv e than those 
probed in current surveys (Barkana & Loeb 2006). 

One important implication of our best-fit duty cycle is that 
~80% of dark matter halos of a given mass are not traced 
by LBGs. The "missing" dark matter halos may have gone 
through bursts of star formation at earlier times and may be 
currently quiescent. However, this does not mean that 80% 
of the stellar mass is missing from observations at z — 6. 
Rather, the gas in the "missing" dark matter halos may not 
have cooled sufficiently to be forming stars rapidly enough to 
be selected as LBGs and thus may not be significant repos- 
itories of stellar mass. The first option suggests that there 
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FIG. 5. — The LBG luminosity function at z=6 (solid line) and z=7.6 (dotted 
line) obtained using the model parameters that maximize the likelihood at 
z=6. The solid and open ci rcles corres pond to observed LBG abund ances 
at z=6 (Bouwens et al. 2006) and z=7.6 I Bouwens & Illingworth 2006). The 
two datapoints at z=7.6 (offset horizontally for clarity) correspond to mor e 
and less conservative selections of z-drops in Bouwens & Illingworth 12006). 
While there is not yet much data at z=7.6, the existing data at these two 
redshifts can be fit without any evolution in the fit parameters and erx> 

may have been a significant amount of star formation at ear- 
lier times. This is evidenced by observations of LBGs at z — 6 
with stellar ma sses as great as a few x 10 10 M Q and ages of 
200-700 Myr dEvles et alj|200ll2006l) . Given the current ob- 
served star formation rate of these galaxies, the past star for- 
mation rate had to have been higher at earlier times. Taken 
together, these observations and the 20% duty cycle inferred 
from the luminosity function of the LBGs suggest that the pur- 
ported deficity o f ionizing photons compared to taht requir ed 
for reionization dBunker et al.| [2004b; Bo uwens et al. 2006) at 
Z ~ 6 could be accounted for by earlier star formation. 

4.2. Lya emitters 

We now use a procedure similar to that described for the 
LBGs to model the LAEs. The key difference is that we must 
also consider the fraction of Lya photons that escape from the 
galaxy and IGM, T a . We generate a grid of models at z = 5.7 
and z = 6.5 with the duty cycle, eoc ranging from 10~ 3 and 1 
and the product of the star formation efficiency and Lya es- 
cape fraction, LT Q spanning between 10~ 3 and 1. Each model 
is compared to the observed abundances, and the likelihood 
is then determined for each model in an identical fashion to 
that discussed for the LBGs. We first perform this procedure 
for our simple model and then examine it in the context of a 
model including supernova feedback. 

Likelihood contours are presented in the top panel of 
Figure [6] for z = 5.7 (solid contours) and z = 6.6 (dot- 
ted contours). As with the fit to the LBGs, there ex- 
ists a strong degeneracy between epc and f+T a . The best- 
fitting parameters (with associated one-sigma uncertainties) 
are (e DC ,f*T a ) = (0.0016^8^,0.0056^^) at z = 5.7 and 
(e D c,f*TJ = (1.0!™ 0.063™) a t z = 6.5. The z = 5.7 data 
are significantly better fit (factor of 4 greater maximum like- 
lihood at z = 5.7) by the advanced model including supernova 
feedback (bottom panel of Figure|6]l. The best-fitting parame- 
ters at z = 6.5 remain unchanged, while those at z = 5.7 change 
slightly: (e DC ,f*T tt ) = (0.040^,0.016^17). Since the 
model with supernova feedback provides a better fit to the 
z = 5.7 data, we focus our discussion on the model parameters 
derived in this fit, rather than the most simple model, in our 
discussion below. 

The best-fitting luminosity functions are plotted over the 
observed abundances in the top panel of Figure [7] (simple 



0.100 




1.000 



H | 0.010 r 



0.001 1 I 

0.01 0.10 1.00 

FIG. 6. — Top: Confidence intervals on free parameters in analytic fit of 
observed LAE abundances at z=5.7 (solid lines) and z=6.5 (dotted lines). The 
likelihood contours are 64% (blue) and 26% (red) of the peak likelihood. The 
parameters that maximize the likelihood are (eDCifsiarT a )=(0. 0016, 0.0056) 
at z=5.7 and (e£>c,fstarT a )=(l. 0,0.063) at z=6.5. Bottom: Same as Figure|6^ 
with the addition of a simple prescription for supernova feedback. 

model) and in the bottom panel of Figure [7] (advanced model 
with supernova feedback). The error bars in these plots in- 
clude both Poisson and clustering variance. Examining the 
confidence intervals and luminosity functions, it seems that 
the data suggest some evolution in the best-fitting model pa- 
rameters between z=5.7 and z=6.5. It is unlikely that the 
evolution is associated with a change in the neutral fraction 
of the IGM because the parameter that is proportional to the 
transmission of Lya photons, f*TL yQ , increases between z=5.7 
and z=6.5, contrary to what would be expected if the neutral 
fraction increased. Taken at face value, the evolution in the 
model parameters suggests that the star formation efficiency 
and lifetime in LAEs increases between z = 5.7 and z = 6.5. 
However, uncertainties in the observations make these con- 
clusions tentative. The LAE luminosity function adopted in 
this paper is based on a photometric sample of objects se- 
lected with a narrowband filter. While many of the most lu- 
minous photometrically-selected objects have been confirmed 
spectroscopically, at lower luminosities the completeness is 
still low. It is possible that there is significant contamination 
from low-redshift line emitters, and hence that the densities at 
these low luminosities are overestimated. If the error bars at 
low luminosit i es are enlarged to reflect this uncertainty, then 
iDijkstra et alJ ([2006) claim the model parameters are consis- 
tent with no evolution between z = 5.7 and z = 6.5 except in 
the underlying halos. Additional spectroscopic observations 
of the lowest luminosity LAEs are clearly necessary to re- 
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FIG. 7. — 7b/>: The Lya luminosity function at z=5.7 (dotted line) and 
z=6.5 (solid line) obtained using the model parameters that maximize the 
likelihood. Open and s olid circles are obs erved abundances at z = 5.7 from 
IShimasaku et alj 1200611 and at z = 6.5 from lKashikawa et all (2006). Bottom: 
Same as above but for best-fitting parameters to model that now includes a 
simple prescription for supernova feedback. 

solve whether the physical parameters of LAEs evolve be- 
tween z=5.7 and z=6.5. 

We now examine the consistency of our best-fit model pa- 
rameters in the context of additional observations of z > 5 
LAEs. Recently, observations have shown that many lu- 
mi nous LAEs a t z — 5 have relatively low stellar masses. 
In iPirzkal et all d2006l) . LAEs at z ~ 5 with L Lya ~ 10 42 - 
10 43 ergs _I have M ste n al - ~ 10 6 - 10 s M Q . Assuming the ra- 
tio of baryons to total mass follows the unversal value J7b/51 m 
and a star formation efficiency of ~ 10%, the observed stellar 
masses suggest halo masses of 0.06-6 x 1O 9 M . According 
to our best-fit models, the halo masses probed by the Subaru 
observations are signficantly greater: at z = 5.7 halo masses 
range from 10 10 to 10 n M Q , while at z = 6.5, the halo masses 
range between 4 x 10 10 and 4x 10 11 M Q . 

Given the low stellar masses observed for these objects, it 
seems that one of our assumptions must be incorrect. Un- 
der our simple model, the stellar mass is given by Miliar = 
f*[fib/^m]Mhai ; m order to decrease the stellar mass for a 
given halo mass, the star formation efficiency must be lower 
than the 10% value we have assumed above. However, the star 
formation efficiency must also satisfy the relation between 
Lya luminosity and halo mass (equation 2) which is con- 
strained via the observed abundance of LAEs as a function of 
luminosity. Hence, if the star formation efficiency decreases, 
either the Lya transmission factor or the ionizing photon rate 
must increase to satisfy equation [2] The latter could occur 
if the IMF of stars was more top-heavy than the standard 
Salpeter form. While a top-heavy IMF would reduce the stel- 
lar mass predicted by our model, it would also decrease the 
observationally inferred LAE stellar masses. These masses 



have been inferred via population synthesis models assuming 
a Salpeter IMF. Top-heavy IMFs have a lower rest-frame op- 
tical stellar mass-to-light ratios then the Salpeter IMF does; 
hence, for a given luminosity, the inferred mass in stars is 
lower than for a Salpeter IMF. Therefore, the stellar masses 
inferred from observations remain at odds with those pre- 
dicted from our models. Alternatively, if the Lya transmis- 
sion fraction, TL yct , is enhanced to account for the lower star 
formation efficiency, the models achieve much better agree- 
ment with the observations. At z = 5.7, our best-fitting mod- 
els have f*TL yQ — 0.04. If the Lya transmission factor is near 
unity, then the star formation efficiency is roughly 4%. In 
this scenario, the stellar masses predicted by the models are 
~ 10 7 — 10 8 Mq at z = 5.7, in much closer agreement with the 
observations. 

The low star formation efficiency and large Lya transmis- 
sion factor of LAEs can be understood physically. In ad- 
dition to showing that the brightest LAEs have low stellar 
masses, observations have also shown that the most lumi- 
nous LAEs are a young population with ages of a few x 10 6 
years (Pirz kal et al.ll2006t iFinkelstein et al.ll2006l) . compara- 
ble to the lifetime of massive stars before they may explode 
as supernovae. Hence, these galaxies are observed at such 
a young stage that they have not had sufficient time to con- 
vert more than roughly 4% of their baryons to stars. More- 
over, they likely have not had enough time to produce a sig- 
nificant amount of dust. This conjecture is corroborated by 
popula tion synthesis mode ling of the observed SEDs of z — 5 
LAEs dGawiser et al.ll2006l) . Without dust to absorb the reso- 
nantly scattered Lya photons, the fraction of Lya photons es- 
caping the galaxy may increase substantially, explaining the 
very large Lya transmission factor needed to fit the observed 
stellar masses. However, the Lya transmission factor is also 
dependent on intergalactic absorption. The typical flux decre- 
ment encountere d by Lya photons in the intergalactic medium 
dFan et al.l l2006) may be substantially reduced in the vicinity 
of LAEs if they reside in groups of galaxies which signifi- 
cantly ionize their surroundings, allowing the Lya photons to 
escape out o f resonance before encountering neutral hyd rogen 
in the IGM (IWvithe & Loebll200l iFurlanetto et aDl20Q6h . 

4.3. Comparison of LAEs and LBGs 

Recent observations at z — 5 suggest that LAEs may differ 
from LBGs in the ir typical stellar m ass and ages. Observa- 
tions presented in lEyles e t al. (2006) have shown that z — 6 
LBGs are a composite population of galaxies, some with low 
stellar masses (~ 10 s M Q ) and some with high stellar masses 
(~ 10 1() M Q ). These stellar masses emerge from our model 
given the best-fit star formation efficiency and halo masses 
probed by the observations. While the uncertainties are still 
significant, there appears to be a weak c orrelation betwee n the 
stellar mass and age of z — 6 LBGs (lEyles et al.ll2006l) : the 
most massive galaxies appear to have a significant population 
of old stars (ages up to ~ 700Myr), while the less massive 
objects appear to be much younger (ages of ~ 60 Myr). The 
model considered in this paper fixes the star formation duty 
cycle to be constant with halo mass and hence cannot confirm 
this observational inference. The best-fitting duty cycle for 
LBGs suggests an age of 200 Myr, which is roughly in the 
middle of the range of lifetimes expected of the LBGs. 

LAEs at z — 5 also appear to be a composite population 
spanning a range of masses and ages, but it appears that their 
typical ages and stellar masses are systematically lower than 
LBGs. A correlation exists between the equivalent width 
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EW) of the Lyg line and the galaxy age and stellar mass 
Finkelstein et al.ll2006l) . The highest EW lines exhibit the 
lowest ages (few million years) and stellar masses (10 6 - 
10 7 Mg); these objects are very faint con tinuum emitters 
and h e nce are not selected in LBG surveys ( Fink elstein et al.1 
I2006t iPirzkal et al.1 120061) . LAEs with lower EWs have 
larger infer red ages (40-200 Myr) with stellar ma sses up to 
~ 10 10 M dFinkelstein et al.l2006tlLai et al.l2007l) ; however, 
these ages and masses do not reach the values as large as those 
seen in LBG surveys at these redshifts. Finally, the SEDs of 
the young LAEs at z — 5 sug gest that there is little exinction 
from dust in these galaxies (IPirzkal et al.l200a) . As explained 
in the §4.2, in order for our model to fit the observations de- 
scribed above, the star formation efficiency of the high EW 
LAEs must be low (due to their extreme youth) and transmis- 
sion of Lya photons through the host galaxy and IGM must 
be very high. 

It is intriguing to consider the fate of the high EW LAEs. In 
one possible scenario, star formation continues, depleting the 
gas content and increasing the total stellar mass of the galaxy. 
The dust content of the galaxy begins to increase. The dust 
absorbs the resonantly-scattered Lya photons, although the 
Lyq EW might be enhanced depen ding on geometric details 
dNeufeldll99 UlHansen & 0120 06). If the star formation rate 
remains high, these objects could continue to be observed as 
LBGs. As more of the gas is converted to stars, the star forma- 
tion efficiency will increase. Within our model we indeed in- 
ferred that the average star formation efficiencies of the z — 6 
LBGs is on order 10%. 

Alternatively, the high EW LAEs could eject their gas via 
feedback from supernovae explosions or quasar activity. As 
soon as the gas density is significantly diluted, the recombina- 
tion rate decreases and there is little emission of Lya photons. 
The galaxies would continue to be selected as LBGs as long 
as massive stars remain; however, without gas star formation 
eventually ceases, leaving the galaxies quiescent and with a 
low stellar mass. Such objects would no longer be detected as 
either LBGs or LAEs. 

5. INTERPRETATIONS OF OBSERVATIONS AT Z ~ 7 - 10 

We now use our model fits, noting the uncertainties, to 
make predictions for both LBGs and LAEs observed at z >7. 
We will assume no evolution in the model parameters to deter- 
mine what redshift trends are expected solely from the natural 
growth of dark matter halos over the era 5 < z < 10. A mod- 
est amount of observational data is available for the z > 7 uni- 
verse and we will address this to see if it is consistent with no 
evolution. Necessarily this discussion will be tentative given 
the considerable uncertainty about the validity of the observa- 
tions. Most of the sources claimed to lie beyond z — 7 have 
no convincing spectroscopic identification. 

5.1. Lyman-Break Galaxies 

Preliminary constraints are no w available on the 
abundance of LBGs at z > 7 (jBouwens et al] 120051: 
iBouwens & Illingwo rth 2006; iRichard et all 120061) . LBGs 
selected as z-band dropouts are considered to have a 
mean redshift z — 7.4. The most recent compilation from 
fields with deep HST-NICMOS data (e.g. GOODS, UDF, 
UDF-P) includes one candidate in the most conservative 
selection and 4 candidates in a more aggressive selection 
(IBouwens & Illingworthl 120061) . After comparing the ob- 
served abjmd^nc£_of_£andjdate 7.4 LBGs to those at 
Z — 6, IBouwens & Il lingworth d2006l) suggest that there is a 
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FIG. 8. — Comparison of model z = 9 LBG luminosity function with 
constraints from observations. The observed abunda nce of LBGs presented 
in the cluster lensing survey of Richard et al. (2006), denoted by solid cir- 
cles, is greater than the uppe r limit found in the Hubble Ultra Deep Field 
(open circle), as pre sented in|Bouwens et al. 12006). The error bars on the 
Richard et al. 1 2003) data are large but may overestimate the true uncertainty 
(see §5.1). The luminosity function obtained by assuming the duty cycle and 
star formation efficiency remain constant between z = 6 and z = 9 (dotte d 
line) is underpredicts the large abundances observed in Richard et al, 120061) . 
Reconciling these data points with our models without resorting to more top- 
heavy IMFs require a star formation efficiency near 100% (solid line). 

rapid assembly of the most luminous star forming systems 
between z — 6 and z — 7.4. 

By extrapolating our star-formation model to z — 7.4 and 
holding the parameters fixed at their best-fit z — 6 values, we 
can determine whether the claimed rapid assembly of lumi- 
nous galaxies requires evolution in the star formation effi- 
ciency o r duty cycle. If all four candi date z-dropouts iden- 
tified in IBouwens & Illingworthl (120061) are at z — 7.4, then 
the observed evolution in the abundance of luminous galaxies 
can be explained simply by evolution of the host dark mat- 
ter haloes (Figure |4j. However, if only one of the candidate 
z-dropouts is at z — 7.4, then our simple model does permit 
some evolution in either the star formation efficiency or duty 
cycle in the 200 million years between z — 7.4 and z — 6. 
Such evolution could be triggered, for example, by the photo- 
ionization heating of the intergalactic medium at the end of 
reionization and the corresponding change in its accretion rate 
onto galaxies. 

At z — 10, several candidate LBGs (sele cted as J-band 
dropo uts) have been identified in the UDF (IBouwens et all 
120051) and in cluster lensing fields (IRichard et al.ll2006h . With 
the addition of deep er optical data, two of the three candidate 
z ~ 10 LBGs from Bouwens et al. (2005) are now known to 
be at lower-redshift (R.J. Bouwens 2006, private communica- 
tion). In this case, the abundance implied if the remaining one 
candidate is at z — 10 is consistent with hierarchical growth 
(Figure |8). 

In contrast, th e abundance o f less luminous z — 9 candi- 
dates located in IRichard et alJ (120061) is significantly larger 
than expected (Figure [8}. However, if clustering fluctuations 
are included in the uncertainties, then the abundances de- 
riv ed are formally consi stent with the lower density implied 
by Bouwens et al. (2005). We note that the 1-sigma error bars 
on the Ric hard et al.l d2006l) data include a large contribution 
(typically a factor of 3) from the uncertainty in the complete- 
ness correction. This uncertainty is difficult to determine ac- 
curately; if overesti mated by as lit t le as ~ 20%, the one-sigma 
uncertainties on the IRichard et al.l (|2006l) datapoint s would no 
longer be consistent with the Bou wens et al.l d2005l) observa- 
tions. More clusters must be studied to verify the large density 
of lensed z — 9 candidate LBGs. If the large abundances are 
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representative, the iRichard et al.l ((2006) observations require 
either significant evolution in the parameters of the simple 
model we have adopted or a top-heavy stellar IMF. Holding 
the duty cycle fixed at its z — 6 value, an implausible star for- 
mation efficiency of 100% is required to explain the observed 
abundance (Figure H). 

The emerging physical picture describing the evolution of 
LBGs at z — 7- 10 is still somewhat tentative. Nonetheless, 
results from the UDF suggest that evolution in the abundance 
of luminous LBGs in the 500 Myr between z — 6-10 can be 
largely explained by the hierarchical assembly, i.e. without 
any evolution in the star formation efficiency or duty cycle. 
In contrast, the high abundance of l e ss lum inous z — 10 can- 
didates suggested by IRichard et alj d2006l) may require sig- 
nificant evolution in either the stellar IMF or the star forma- 
tion efficiency. Sim ilar conclus i ons ar e reached using semi- 
analytic models in Sam ui et ail d2006l) . Additional observa- 
tions are required to confirm the large density observed in 
IRichard et al.l ((2006) is robust and to reconcile these poten- 
tially differing pictures. In §6, we will use our model to pre- 
dict the ability of future surveys to detect starburst galaxies at 
z~7-10. 



5.2. Lya Emitters 

The first results are also now available from Lya surveys at 
z ~ 9. As with the LBGs, the high redshift results are seem- 
ingly contradictory. Willis et al. (2005) and Cuby et al. (2006) 
find no Lya emitters in narrowband surveys centered at z = 8.8 
with the Very Large Telescope (VLT). However, both surveys 
are only sensitive to the brightest LAEs (> 3 x 10 42 erg s" 1 in 
I Willis & Courbinl2005l and >10 43 erg s" 1 in lCubv et alj2006h 
over modest comoving volumes (870 Mpc 3 and 5000 Mpc 3 , 
re spectively). 

IStark et all (120071) conducted a cluster lensing survey for 
LAEs at z = 8.5 - 10.4. The magnification provided by the 
clusters allows significantly less luminous LAEs to be de- 
tected (>10 41 erg s" 1 ), albeit over a much smaller comoving 
volume (~ 30 Mpc 3 ). Six candidate LAEs were identified 
with unlensed luminosities spanning 2-50xl0 41 erg s _1 ; at 
least two of the six candidates are considered likely to be at 
z ~ 9 following additional spectroscopy which casts doubt on 
alternative low-redshift explanations for the J-band emission 
features. 

In Figure|9] we use our model to compute luminosity func- 
tions of LAEs at z — 9 assuming the duty cycle, star formation 
efficiency, and Lya escape fraction remain fixed at either their 
best-fit z = 5.7 or z = 6.5 values. As the narrowband observa- 
tions refer to a single redshift, we use the mass function at 
z = 8.8 to compute the lumin osity funct i on. F or the cluster 
lensing survey presented by IStark et all (120071) . as the halo 
mass function evolves significantly over the redshift range 
sampled, we compute the average halo mass function between 
z = 8.5 and z = 10.4, weighting mass by the relevant sensitiv- 
ity function. We generated luminosity functions using both 
mass functions. While the resulting luminosity functions are 
marginally different, the net results described do not change; 
hence, for the sake of clarity, in Figure [9] we only overlay the 
luminosity function from the z = 8.8 mass function on the data 
points from the three surveys described above. 

The results suggest that, for luminous LAEs, current sur- 
veys do not yet have the combined sensitivity and depth to 
detect an y sources at z — 9 . Alth oug h the upper l imits pre- 
sented in lWillis & Courbinl (120051) and lCubv etail (120061) are 
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FIG. 9. — Comparison of model z — 9 LAE luminosity functio n with con- 
straints from observations. Top: If all the z — 9 LAE candidates in l Stark et all 
1 2007) are at high-redshift, then the best-fit model parameters at z = 5.7 (dot- 
ted line) and z = 6.5 (solid line) are inconsistent with their values at z = 9, 
Bottom: Model luminosity functions assuming an extreme top-heavy IMF 
characteristic of Pop-Ill stars with model parameters fixed at their z = 6.5 
(dotted line) and a star formation efficiency of 100% (with duty cycle fixed at 
its best-fit z=6.5 value, solid line) overlaid upon observational constraints. 

consistent with our expectations, those surveys only rule out 
the possibility that the density of luminous LAEs decreases 
in the time interval between z — 9 and z — 6. On the other 
hand, for l e ss lum inous LAEs, if all six of the candidates in 
IStark et ail d2007l) are at high-redshift and the inferred abun- 
dances are representative, significant evolution is implied in 
the model parameters. As with the lensed LBGs at z — 9, ei- 
ther a high (~ 100%) star formation efficiency (for a fixed 
duty cycle) or a top-heavy IMF of stars (Figure |9]l would be 
required. An important caveat is the uncertainty caused by 
cosmic variance, which we do not include in the error bars 
in Figure|9] Even for a significantly more ambitious spectro- 
scopic lensing survey, the fluctuations expected from large- 
scale structure are > 100% (Figure |3); hence, it is possible 
that the candidate LAEs discovered in lStark et al.l d2007l) may 
trace an extreme overdensity in the underlying mass distribu- 
tion, in which case the derived densities may be larger than 
the cosmic average at that epoch. Clearly, more clusters must 
be observed. If only the two prime LAE candidates are at 
high redshift 3 , the derived abundances are formally consistent 
with the predicted z — 9 luminosity function due to the large 
Poisson fluctuations. 

6. IMPLICATIONS FOR FUTURE SURVEYS 



3 The fact that the Poisson errors are 100% for the two datapoints in Figure 
^corresponding to the case that two candidates from Stark et al. 1 2007) are at 
z ~ 10 may seem counterintuitive given that there are two sources; however 
this arises from strongly-varying limiting (detectable) source lum inosity over 
the fie ld of view due to the cluster magnification. For details, see S tark et aT] 
<200l) . 
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Over the next few years, several new instruments will be 
placed on current ground-based telescopes and Hubble Space 
Telescope motivated in part by extending the search for z >7 
sources. These surveys offer the exciting possibility of char- 
acterizing the assembly of the first galactic sources through 
the era of reionization. With our model, calibrated by the data 
available at z — 6, we can consider the optimal volume and 
depth necessary for detecting star-forming sources of various 
luminosities to z — 10-20. Thus it is hoped our model can 
assist in guiding the design of future instruments and surveys. 

In the following subsections, we will evaluate several 
benchmark imaging and spectroscopic surveys. We consider 
a large, blank-field, narrowband survey for LAEs at z — 8 and 
Z— 10 in §6.1, and discuss two different surveys for dropouts 
at z — 7.5 and z — 10 in §6.2. Finally, in §6.3, we examine 
the efficiency of a variety of lensing surveys for LAEs and 
LBGs. Since some of these proposed surveys reach to signif- 
icantly lower luminosities than the current surveys discussed 
in §4 and §5, we examine the effects of feedback on the pre- 
dicted counts, as well as the gains of adaptive optics given our 
predicted sizes for faint LAEs. 

6.1. The Dark Ages z Lyman-alpha Exlorer: LAEs at 
z-7-10 

The Dark Ages z Lyman-alpha Explorer (DAzLE) is a nar- 
rowband imager on the VLT which aim s to detect Lya emit- 
ters at 6.5 < z < 12 dHorton et al.ll2004l) . DAzLE has recently 
observed two pointings of GOODS-South in two filters corre- 
sponding to Lya redshfits of ~ 7.7 and ~ 8.0 (R. McMahon 
2006, private communication). The observing sequence for 
DAzLE involves alternating between two narrowband filters 
with slightly different central wavelengths. A composite im- 
age is made of all of the subexposures in each filter, resulting 
in two "subsurveys" slightly offset in redshift space. Sub- 
tracting the two composite images removes continuum soures, 
thereby allowing candidate LAEs to be identified. 

In ten hours of integration, DAzLE is expected to reach a 3a 
sensitivity of 2 x 10 " 18 erg cm" 2 s" 1 in the differenced image 
dHorton et al.ll2004b . At z = 7.7, this corresponds to a limiting 
LAE luminosity of 1.5 x 10 42 erg s . As a benchmark sur- 
vey with this instrument, we consider a four position mosaic 
(i.e. 4x 6/83x6/83). The total comoving volume sampled 
in four pointings of the two filters at z = 7.7 is ~ 6900 Mpc 3 . 
A simple extrapolation of our model suggests that a comov- 
ing volume of 1100 Mpc 3 is necessary to detect one LAE at 
z = 7.7. Thus, in this proposed survey, 6-7 LAEs would be 
detected with DAzLE assuming no rapid evolution in the star 
formation efficiency, duty cycle or IGM transmission between 
z ~ 6.5 and z — 7.7. 

However, although a promising survey in terms of likely 
detections, the uncertainties are considerable. Clustering fluc- 
tuations are 40-50% (for the best-fit z=5.7 and 6.5 model 
parameters, respectively) at the limiting luminosity of 1 .5 x 
10 42 erg s" 1 . The Poisson fluctuations expected for a 7200 
Mpc 3 survey are similar (~ 40%). It is prudent to con- 
sider what effects such large fluctuations would have on at- 
tempts at using the DAzLE results to constrain the progress of 
reionization via the evolution of the LAE luminosity f unction 
(iMalhotra & PJioadsll2lM 120061: iDiikstra et^l2006h . While 
our model predicts that 6-7 sources should be detected assum- 
ing only evolution in the underlying halos, the large (~ 60%) 
expected field-to-field variations imply that the source counts 
could vary significantly from the predicted value. Ignoring 
additional complications on the tranmission of Lya photons 
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FIG. 10. — Top: Extrapolation of Lya luminosity function to z — 1-1 (dotted 
red curve) and 8.3 (dashed purple curve) fixing star formation model param- 
eters at their best-fit z. = 6.5 values (§4). The predicted luminosity function 
suggests that DAzLE should detect 6-7 z — 7.7 sources in a four pointing 
mosaic with total integration time of ~ 40 hours. A lensing survey utitlizing 
F2T2 could detect up to 8 z — 8.3 sources in ~ 200 hours of observations. 
Bottom: Extrapolation of Lya luminosity function to z — 10 fixing star for- 
mation model parameters at their best-fit z = 6.5 values. Almost 200 hours 
are required to detect a single z, ~ 10 source with DAzLE under these as- 
sumptions. The estimated performance of lensing surveys with a multi-object 
spectrograph such as MOSFIRE or a tunable narrowband filter (F2T2) sug- 
gest that lensing is an efficient means of detecting LAEs at z. — 7 — 10, but 
supernova feedback could drastically reduce the number of LAEs detected. 

through the IGM from galax y groups (IWvithe & Loebl 120051 
iFurlanetto et al.ll2004l 120061) and peculiar velocities (Dijkstra 
& Loeb 2006, in preparation), if between two and ten LAEs 
are detected with DAzLE in GOODS-S, little information can 
be reliably deduced on the evolution of IGM. Holding the 
duty cycle and star formation efficiency fixed, it would re- 
quire a >60% decrease in the Lya transmission factor, T a , 
for only one z = 7.7 LAE to be detected towards GOODS-S. 
Hence, these results suggest that, with currently feasible sur- 
vey geometries and instruments, this method of constraining 
reionization will only be effective if the IGM evolves rapidly 
(>60%) over a short redshift interval. 

DAzLE is able to search for LAEs out to z — 10. There 
is a large gap in the atmospheric OH forest between 1.3325 
and 1.3401 ^m, corresponding to a LAE redshift win- 
dow of Az=0.06 centered at z = 9.99. A similar window is 
located at z = 9.91. A single DAzLE pointing at this redshift 
(assuming a filter width of ~ 10A) samples a comoving vol- 
ume of ~ 1250 Mpc 3 summed over both filters. Extrapolating 
the model luminosity function to this redshift with the param- 
eters fixed at their best-fit z = 5.7 values, it would take just 
over a volume of 2500 Mpc 3 (two pointings) to detect one 
LAE brighter than 10 42 erg s" 1 . At z = 9.99, this luminosity 
limit corresponds to a limiting flux of 8 x 10~ 19 erg cm" 2 s" 1 . 
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Reaching this sensitivity would require 70-80 hours for each 
pointing (adjusting for the expected atmospheric transparency 
in the wavelength interval of the observations); hence almost 
200 hours would be required to detect a single z — 10 source. 
Altern atively, if the LAE candidates identified by Sta rk et all 
(2007) are at high-redshift, then > 100 LAEs brighter than 
10 42 erg s" 1 would be expected in the single 70-80 hour DA- 
zLE pointing. 

In summary, DAzLE may well detect many z=7.7 sources 
but spectroscopic confirmation will be a challenge. However, 
even in the most ambitious surveys we can currently contem- 
plate, the contribution of these sources to reionization will 
be seriously limited by the expected clustering fluctuations. 
Surveys at higher redshift will be much more demanding. 
Few sources are expected within reasonable exposure times 
at z — 9-10 unless there is significant evolution in the LAE 
model parameters between 6< z < 10, as may be possible only 
if all of the fainter lensed LAEs detected bv lStark et all (120071) 
are at high redshift. 

6.2. Imaging Surveys for LBGs at z — 7 — 10 

Selection of z > 7 LBGs will be greatly aided by the 
new generation of large-format near-infrared detectors. The 
Subaru M ulti-Object Infrared Camera and Spectrograph 
(MOIRCS, llchikawa et alj 120061) offers imaging and spec- 
troscopic capabilities over a large 4' x 7' field of view. In 
2008, the Wide Field Camera 3 (WFC3) is scheduled to be 
installed on HST; this near-IR camera will offer imaging in 
a 127" x 137" field of view. We thus consider likely sur- 
veys with MOIRCS and WFC3. First we consider a MOIRCS 
wide-field mosaic whose sensitivity is arranged to match that 
of the near-infrared observations of GOODS-S, but with an 
area twice as large. We then consider a single, ultra-deep 
pointing with WFC3. 

To observe an area twice as large of GOODS-S would 
require eleven MOIRCS pointings. Selecting reliable z- 
dropouts (z — 7.5) and J-dropouts (z — 10) requires deep z, 
J, H, and K-band data. We assume the near-infrared data is 
roughly similar in depth to GOODS, with 5-er point sources 
sensitivies of z'=26.6, J=25.8, H=24, and K=24. These lim- 
its are optimized to the selection cr iteria of z'-drops (z' - 
J >0.8.lBouwens & Illingworfhl2006h and J-drops (J-H>1.8, 
iBouweiis et al.ll2005h . The K-band limit is chosen with the 
goal of detecting J-drops in a second filter to help remove 
false-positives. The limiting J- and H-band limits correpond 
to limiting star formation rates of 18 and 146 M yr _1 at z — 
7.5 and 10, respectively. Reaching J=25.8 with MOIRCS 
requires ~ 12 hours of integration assuming 0"5 seeing (J. 
Richard 2006, private communication). Significantly less 
time is needed to reach the desired sensitivies in H (1.8 hours) 
and K (40 minutes). The z'-band observations would be most 
efficiently performed with the Suprime-Cam on Subaru. The 
entire area could be covered in one pointing with Suprime- 
Cam, requiring 4.7 hours of integration. In total, 163 hours 
would be required for the observing sequence. We consider 
this a practical, albeit ambitious, program. 

Each MOIRCS field of view samples a comoving area of 
10.6 Mpc x 18.6 Mpc at z ~ 7.5 and 11.4 Mpc x 20.0 Mpc 
at z — 10. We assume the redshift probability distribution of 
z-drops and J-drops is a Gaussian with a standard deviation of 
(j;=0.5 and a mean of /i z =7.5 and 10 for z-drops and J-drops, 
respectively. We normalize the probability distribution so that 
the maximum completeness is 50% at the mean redshift of the 
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FIG. 1 1. — Predicted LBG luminosity function at z — 7.5 assuming star 
formation model parameters are fixed at their best-fit z — 6 values (§4). The 
absolute magnitude at 1350A is plotted horizontally along the bottom of the 
plot, and the corresponding J-band apparent magnitude is plotted along the 
top of the plot, assuming z = 7.5. The efficiency of several mock surveys for 
Z ~ 7.5 sources is overplotted. 

survey. While the completeness is certainly higher for objects 
that are much brighter than the sensitivity limit of the survey, 
monte carlo simulations suggest that 50% is a reasonable av- 
erage fo r entire population which is dominated by the faintest 
objects (IStark et al.ll2006l) . We compute the effective radial 
distance sampled by the surveys by integrating the normal- 
ized redshift probability distribution over a distance spanning 
Az=2 in redshift and centered at /i z . Over eleven pointings, 
this corresponds to a total comoving volume of 4.4 x 10 5 Mpc 3 
at z ~ 7.5 and 3.5 x 10 5 Mpc 3 at z ~ 10. 

In Figures QT| and [T2] we extrapolate the best-fit luminos- 
ity function from z — 6 to z — 7.5 and z — 10 allowing only 
for evolution in the dark matter mass function. The comoving 
number density of detected LBGs at the sensitivity limit of 
the MOIRCS survey is 1 .2 x 10" 5 Mpc" 3 mag" 1 at z ~ 7.5 and 
1.3 x 10~ 10 Mpc" 3 mag" 1 at z — 10. Integrating over the en- 
tire magnitude range, 1-2 sources would be detected brighter 
than J=25.8 at z — 7.5 and 4x 10~ 6 sources would be detected 
brighter than H=24 at z — 10. Clearly, neither the z — 7.5 
or z — 10 MOIRCS mock-survey is very efficient at detecting 
high-redshift sources. Even with the relatively large areal cov- 
erage provided by MOIRCS, the integrated near-infrared sky 
is simply too bright in broadband filters to reach the sensitivity 
limits necessary to detect an abundant population of z > 7- 10 
LBGs; with current technology, conventional ground-based 
surveys are much better off searching for LAEs using nar- 
rower filters tuned in between the bright sky lines. 

As with the MOIRCS observations, the mock WFC3 ob- 
servations require deep z, J, and H-band coverage to select 
z- and J-band dropouts. WFC3 offers a significant improve- 
ment in both throughput and areal coverage with respect to 
NICMOS, so the camera is ideal for selecting J-drops. By 
contrast, WFC3 is not more efficient than ACS in the z-band. 
Therefore, WFC3 is not particularly well-suited for conduct- 
ing a new z-drop survey that is significantly deeper than UDF 
We thus consider its efficiency at detecting J-dropouts. 

We adopt limiting magnitudes that are ^ 1 and ~2 magni- 
tudes deeper than the UDF in Hi6o and Jno, respectively, cor- 
responding to 5(7 sensitivities of Jhq=31.3 and Hi6o=29.5 for 
point sources (assuming an 0."4 diameter aperture). With a 
J-band 5cr sensitivity of 31.3, J-drops can be s elected as faint 
as Hi6 o=29.5, using the selection criteria of iBouwens et alJ 
(2005). The estimated star formation rate for z — 10 LBGs 
with Hi6o=29.5 is 0.9 M Q yr _1 . Based on the anticipated 
WFC3 performance, reaching such sensitivities in a single 



13 



m 1350 

21 22 23 24 25 26 27 28 29 30 31 32 



c 



10-' 





WFC3 LENS ■ 










MOIRCSi 




1 









-26 -24 -22 -20 -18 -16 

M 1350 



FIG. 12. — Predicted LBG luminosity function at z — 10 assuming star 
formation model parameters are fixed at their best-fit z — 6 values (§4). The 
absolute magnitude at 1350 A is plotted horizontally along the bottom of the 
plot, and the corresponding H-band apparent magnitude is plotted along the 
top of the plot, assuming z = 10. The efficiency of several mock surveys for 
Z ~ 10 sources is overplotted. 

pointing would take 301 hours in Jno and 48 hours in Hi6o 
resulting in 349 hours of total integration 4 . Since K-band 
observations are not practical with WFC3, we rely only on 
an H-band detection for J-drop selection. Each WFC3 field of 
view covers a comoving area of 6.0 Mpc x 6.5 Mpc at z — 10. 
Assuming a redshift probability distribution identical to that 
described above for z-drops and J-drops with MOIRCS, we 
find that the single WFC3 pointing will sample a comoving 
volume of 5410 Mpc 3 at z ~ 10. 

The predicted comoving number density of detected J- 
drop LBGs at the sensitivity limit of the single deep WFC3 
pointing is 7.2 x 10~ 4 Mpc" 3 mag" 1 at z — 10, assuming no- 
evolution in the star formation efficiency and duty cycle. In- 
tegrating over the entire magnitude range, 1 -2 sources would 
be detected at z — 10. 

In summary, existing ground-based imagers are not well- 
equipped for detecting LBGs at z — 7 — 10, largely because 
the integrated near-IR sky is too bright to reach the sensitiv- 
ity limits necessary to detect sources at z > 7 in a reasonable 
amount of time. This problem could be lessened if the field 
of view was much bigger, allowing the the most luminous and 
rare LBGs to be detected. Even from space (using WFC3 on 
HST), detecting LBGs at z — 7- 10 will not be trivial, requir- 
ing hundreds of hours to detect a single LBG at z — 10. 

6.3. Lensing Surveys for Star-forming Galaxies at z — 7 — 10 

Currently, strong lensing surveys for high-redshift galax- 
ies could be contemplated for about 20 galaxy clusters for 
which there are well-defined mass models; these mass mod- 
els are essential in accurately defining the spatial distribution 
of magnification. However, ongoing HST surveys will signif- 
icantly increase the number of suitable galaxy clusters (e.g. 
Ebeli ng et al.ll2003h . Here, we examine the efficiency of sur- 
veys for lensed LBGs and LAEs assuming a larger sample of 
clusters and more efficient instruments that will soon become 
available. We consider both a n extension of the longslit spec- 
troscopic survey discussed in (IStark et al.l2007l) and as well as 
an imaging campaign to identify lensed z-drops and J-drops. 

In the next several years, a number of near-IR multi-object 

4 E stimated sensitivie s for WFC3 are listed at 

http : / /www . stsci . edu/hst/wf c3 /documents /handbook/ 
cyclel6/wf c3_cycl66 .html 



spectrometers will be installed on 8-10 meter class telescopes, 
offering significant gains in sensitivity and field-of-view over 
current near-IR instruments. One such instrument is the 
Multi-Object Spectrograph For Infra-Red Exploration (MOS- 
FIRE, PI: I. McLean & C. Steidel) on Keck I. MOSFIRE will 
utilize a configurable slit unit allowing up to 45 slits, each 7. "3 
in length, within the 6/14x6/14 field of view. While tilted 
and curved slits are not possible with MOIRCS, the config- 
urable slit unit will still allow for more areal coverage of high 
magnification regions than with NIRSPEC. In addition, MOS- 
FIRE will be more sensitive than NIRSPEC, reaching 5cr line 
flux sensitivities of 1 x 10~ 18 erg cm" 2 s" 1 in 4 hours, assuming 
an unresolved line with R=3270 5 . Given the specifications of 
MOSFIRE, it is optimal to concentrate observations to clus- 
ters with vary large and well-determined critical lines (e.g. 
Abell 1689, Abell 1703, Abell 2218). While the exact magni- 
fication distribution provided to background sources may very 
slightly from cluster to cluster, it is reasonable to expect 80% 
of the MOSFIRE survey area to be at a magnification of 20, 
15% of the area to have a magnification o f 10, and 5% of th e 
area to be magnified by only a factor of 5 dStark et al.ll2007l) . 

We consider a MOSFIRE spectroscopic lensing survey for 
LAEs both at z = 7.0-8.3 (Y-band) and at z = 8.5 - 10.3 
(J-band). Assuming ~ 15% of the line-of-sight distance is 
lost due to bright OH lines, this corresponds to comoving 
radial distances of 353 and 334 Mpc, respectively. Assum- 
ing 10 clusters are observed for 8 hours each in the Y-band, 
the resulting survey volume is ~ 73 Mpc 3 for LAEs brighter 
than 10 40 5 ergs"'. For the mock J-band survey, we assume 
longer integration times (20 hours) over eight clusters, giv- 
ing an identical integration time to the mock DAzLE z = 9.9 
survey considered previously. This results in the coverage of 
~ 58 Mpc 3 for z ~ 10 LAEs brighter than lO^ergs" 1 . Ex- 
trapolating our LAE model to z = 7.5 and to z = 10 holding the 
model parameters fixed at their z = 5.7 values, we predict that 
the MOSFIRE survey should detect 9 sources brighter than 
10 405 ergs _1 at z = 7.5 and 4 sources brighter than this limit at 
z ~ 10. If supernova feedback decreases the star formation ef- 
ficiency in a manner described in §2, then the predicted num- 
ber of LAEs would drop drastically. The z — 7.5 survey would 
potentially detect 1-2 LAE brighter than 10 40 ' 5 ergs _1 while the 
z ~ 10 survey would not detect any sources. Hence, the suc- 
cess of the mock lensing survey depends strongly on whether 
supernova feedback decreases the efficiency of star formation 
in low-mass dark matter halos. Given the possible efficiency 
with which the lensing survey could detect z — 7 - 10 galaxies, 
it is of the utmost importance to observationally constrain the 
effects of supernova feedback on the LAE luminosity function 
at lower redshifts (z — 3 - 6). 

Gravitational lensing can also be very valuable for imag- 
ing surveys for LBGs that are intrinsically fainter than those 
detected in conventional deep surveys (e.g. GOODS, UDF). 
One such survey is currently being co nducted toward six 
galaxy clusters with NICMOS on HST (IStark & Ellisll2006l 
Richard et al. 2007, in preparation). With WFC3, such a 
survey could potentially be conducted much more efficiently, 
allowing many more clusters to be observed. We examine the 
feasibility of a hypothetical WFC3 lensing survey of galaxy 
clusters for z- and J-dropout galaxies. The primary benefit 

5 Sensitivities and additional specifications of the spec- 
trograph are provided in the Preliminary Design Report 
located at |http : / /www . astro . ucla .edu~irlab/| 

most ire/MOSFIRE%20PDR%20Report%20v4 .pdf 
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of a WFC3 survey for z > 7 LBGs is the added throughput 
and field of view in the Jno and Hi6o-bands compared to what 
is available with NICMOS. For each cluster, we assume 5-a 
point-source sensitivities (in an aperture with 0"4 diameter) 
of Z85o=27.4, Jno=28.2, and Hi6o=27.4, requiring ~ 8, 1, and 
1 hour(s) per cluster, respectively. These limits allow z-drops 
to be selected down to Jno=26.6 and J-drops to be selected 
down to Hi6o=26.4 (without considering the effects of lens- 
ing). If we allot 350 hours to this observing program (an iden- 
tical time allocation to the traditional survey discussed in the 
previous subsection), 35 clusters can be observed. The total 
area surveyed would be more than a factor of fifteen greater 
than previous lensing surveys for z > 7.5 LBGs (Richard et 
al. 2007, in preparation). 

Both the survey volume and limiting source luminosity are 
modified by the magnification prov id ed by the foreg round 
galaxy cluster (see lSantos et al.l2004l or lStark et al.l2007l for a 
description). A typical cluster provides a magnification boost 
of x 2,5, 10, and 30 over 92%, 46%, 29%, and 12% of the 
entire WFC3 field of view. Adopting this magnification dis- 
tribution for each of the 20 clusters, at z — 7.5, the survey 
is sensitive to a volume of 2.4x10 s Mpc 3 , 7.6xl0 4 Mpc 3 , 
1.5xl0 4 Mpc 3 , 3.1 xlO 3 Mpc 3 , and 580 Mpc 3 for sources 
brighter than Jno = 26.6, 27.6, 28.6, 29.6, and 30.6 respec- 
tively 6 . At z — 10, the survey probes a comoving volume 
of 1.9xl0 5 Mpc 3 , 6.0xl0 4 Mpc 3 , 1.1 xlO 4 Mpc 3 , 2.4xl0 3 
Mpc 3 and 460 Mpc 3 for sources brighter than Hi6o=26.4, 
27.4, 28.4, 29.4, and 30.4, respectively. Over 35 clusters, 
if the star formation efficiency and duty cycle remain fixed, 
then such a survey should detect 82 z — 7.5 LBGs brighter 
than Ji io=3 1.6 (corresponding to a star formation rate of 0.09 
M Q yr _1 at z — 7.5 for the source assumptions discussed in §2, 
Figure [m> and 6 z ~ 10 LBGs brighter than H 16 o=31.4 (cor- 
responding to a star formation rate of 0.2 Moyr -1 at z — 10, 
Figure [T2l. 

Adaptive optics (AO) provides the possibility of diffraction- 
limited observations from the ground. If the projected size 
of the target objects is small enough, such observations are 
more efficient than non-AO observations because the photo- 
metric aperture can be decreased, thereby allowing signifi- 
cantly less noise for nearly the same amount of flux 7 . One 
planned survey that aims to take advantage of AO is the Gem- 
ini Genesis Survey (GGS ). This survey w ill use a tunable 
Fabry-Perot etalon (F2T2, IScott et al.ll2006l) on Gemini with 
resolution of R=800 to detect lensed LAEs at z ~ 8 - 10. If 
the projected angular size of the sources is less than ~0.03 
arcsec 2 , then GGS should be able to reach a 5a sensitivity of 
3-6 x 10" I8 erg cm" 2 s" 1 in 10 minutes (R. Abraham 2006, 
personal communication). Following the scaling relation de- 
rived in §2, galaxies at z — 10 should have typical sizes of 
0.01 arcsec 2 , small enough to significantly benefit from AO. 
The field of view of F2T2 is 45"x45", ideally suited to imag- 
ing the most highly magnified regions of galaxy clusters. 

We consider a mock GGS survey of 60 clusters; we assume 
each cluster is observed for 5 minutes in 40 different wave- 
length positions between 1.1 /xm and 1.3/xm, allowing the de- 
tection of LAEs at z = 8.1 to z = 9.7. This should take ~ 200 
hours of integration. A typical cluster provides a magnifica- 

6 The quoted magnitudes correspond to the apparent magnitude that a 
source would be observed with if it was not magnified. 

7 The validity of this statement depends on the strehl ratio, which is the ra- 
tio of the peak brightness of the stellar image to that produced by an ideal 
optical system. If the strehl is very low, then the amount of flux in the 
diffraction-limited aperture will be greatly reduced. 



tion gain of x5, 10, and 30 over 80%, 6 9%, and 33% of the 
F2T2 field of view (iRichard et al .1120061) . Taking this as the 
magnification distribution for each of the clusters, the total 
survey volume sensitive to LAEs brighter than 10 41 erg s" 1 
is 30-86 Mpc 3 . The z — 8-8.5 LAE luminosity function (as- 
suming fixed model parameters from z = 5 .7) suggests that 3-7 
sources would be detected at z — 8-8.5. If the effects of su- 
pernova feedback are parameterized as in §2 and §4, the num- 
ber density of low luminosity sources is drastically reduced; 
in this case, no sources would be detec ted in the survey. Al- 
ternatively, if all six candidate LAEs in lStark et alj (120071) are 
at high-redshift, then GGS should detect over 30 z ~ 8- 10 
sources. 

It appears that lensing surveys offer one of the more effi- 
cient means of identifying galaxies at z — 7- 10 since they are 
able to reach sensitivity limits where objects are expected to 
be much more abundant. However, supernova feedback may 
drastically reduce the number of sources detected in these sur- 
veys. Regardless, spectroscopic confirmation of these sources 
will continue to remain challenging until JWST and 20-30 
meter ground-based telescopes become available. 

7. CONCLUSIONS 

We have attempted to empirically calibrate the parameters 
of a simple star formation model using observations of star- 
forming galaxies (both LBGs and LAEs) at z — 6. The error 
budget used in fitting the data takes proper account of both the 
Poisson and clustering variance. We use the calibrated model 
to characterize the physical properties of LBGs and LAEs at 
z ~ 6 and extrapolate it to higher redshifts to make predictions 
for upcoming surveys for galaxies at z — 7- 10. Our primary 
conclusions are as follows: 

1. We have derived accurate formulae for the field-to-field 
variance expected in broadband surveys for LBGs, narrow- 
band surveys for LAEs, and lensing surveys for LAEs. For 
each survey geometry, there exists a cross-over luminos- 
ity L c below which the clustering variance dominates over 
commonly-used Poisson variance. In total, the clustering 
variance accounts for less than 6% error in narrowband sur- 
veys for LAEs in the Subaru Deep Field and 15-20% error 
in the z — 6 surveys for LBGs in UDF The clustering fluctua- 
tions are significantly higher for spectroscopic lensing surveys 
reaching up to 100%. 

2. LBGs at z — 6 are best-fit by a model with a star formation 
efficiency of 13% and a duty cycle of 0.2. The star forma- 
tion efficiency suggests that, on average, 87% of the baryonic 
mass of z — 6 LBGs still remains in the gas-phase. The duty 
cycle indicates that the current burst of star formation has a 
lifetime of 200 Myr, roughly equivalent to the dynamical time 
of virialized halos at z — 6. The duty cycle also implies that 
80% of dark matter halos of a given mass are not traced by 
LBGs. The missing halos could have yet to form many stars 
(due perhaps to inefficient cooling) or could be quiescent after 
experiencing a burst of star formation at earlier times. This re- 
sult suggests that the claimed deficit in ionizing photons from 
luminous LBGs at z — 6 relative to what is required for reion- 
ization may be explained by star formation at higher redshift. 

3. The best-fitting model parameters for LAEs show some 
evidence for evolution between z = 5.7 and z = 6.5. However, 
most likely the evolution is not due to a change in the neutral 
fraction of the IGM since the parameter that is proportional to 
the transmission of Lya photons through the IGM increases 
between z = 5 .7 and z = 6 .5 . Thus, we consider the evolution to 
be tentative because of the large uncertainties in the density of 
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the lowest luminosity LAEs. Additional spectroscopic efforts 
are needed to improve the spectroscopic completeness at low- 
luminosities. 

4. The star formation efficiency of LAEs must be very low 
(~ 1%) in order to reproduce the low stellar masses inferred 
from observations of LAEs at z — 5. Such a low star formation 
efficiency is only possible in the context of the models if the 
Lya transmission factor is near unity. This not only requires 
a large escape fraction of Lya photons from the galaxy, but 
also requires that the photons are not substantially absorbed in 
their path through the IGM. The emerging physical picture is 
that the LAEs with large EWs are young objects (10-20 Myr) 
that have only converted on order 1 % of their baryons to stars 
(and hence remain gas-rich) and have not yet produced much 
dust, allowing a large escape fraction for Lya photons. The 
intergalactic absorption decrement in Lya may be reduced if 
the ionizing luminosity of the LAEs or the galaxy groups in 
which they reside is sufficiently large so as to strongly ionize 
the IGM in their vicinity. 

5. We have attempted to fit preliminary data at z — 7 - 10 by 
extrapolating our model to higher redshifts assuming model 
parameters are fixed at their calibrated z — 6 values. The ob- 
served evolution of LBGs between z — 6 and z — 7.5 can be 
explained largely by changes in the host dark matter halos. At 
z ~ 10 the picture is more complicated. Co nstraints on the 
prese nce of z — 10 LBGs in the Hubble UDF (Bo uwens et alJ 
120051) are explained by the hierarchical growth of dark matter 
halos. In contrast, the larg e abundance of z — 9 LBGs claimed 
in the lensing survey of lRichardet~ai1 d2006) is difficult to 
explain without resorting to a top-heavy IMF or large field- 
to-field fluctuations. The situation is similar for the lensed 
LAEs. Two traditional surveys at z = 8.8 are consistent with 
no-evolution in the model parameters from z = 5.7, while the 



abundances infe rred from candida te LAEs in a spectroscopic 
lensing survey (IStark et aT]|2007l) are only explicable if the 
IMF is top-heavy or if the observations are probing an over- 
density in the underlying mass distribution. 

6. New instruments that will become available on ground- 
based telescopes and HST in the next 2-3 years should greatly 
increase the efficiency with which z — 7 - 8 LAEs and LBGs 
are detected. However, unless there is significant upward evo- 
lution in the luminosity function from z — 6, detecting z — 10 
galaxies via conventional methods will only be feasible if 
heroic efforts are undertaken. With current telescopes, lens- 
ing surveys are potentially better suited to detecting z — 10 
sources depending on supernova feedback. JWST and thirty- 
meter class ground-based telescopes are most likely necessary 
to detect a substantial population of objects at z — 10. 

7. Constraining reionization at z — 6 - 7 via the evolution in 
the LAE luminosity function as probed by future narrowband 
surveys will be complicated by clustering and Poisson vari- 
ance. Given the large fluctuations (> 60% for feasible survey 
geometries and sensitivities) and small number of sources ex- 
pected to be detected (~ 7 assuming no evolution from z=5.7), 
this technique will only be effective if the neutral fraction in 
the IGM evolves very rapidly in the redshift interval between 
z = 6.5 and z = 7.7. 
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TABLE 1 

Variance in High-z Galaxy Surveys 



Field 


Type 


Field of View (arcmin 2 ) 


z 


Llt .c (ergs ') 


Li5oo (erg s 1 Hz ') 




GOODS 


dropout 


160 


6 


N/A 


10 28.0 


0.002 


UDF 


dropout 


11 


6 


N/A 


10 28.2 


0.03 


SDF 


narrowband 


1295 


5.7 


10 41.7 


N/A 


0.001 


Cluster 


spec, lensing 


0.13 


8.5-10.4 


N/A 


N/A 


N/A 



NOTE. — The field of view for the spectroscopic lensing survey is in the source plane, assuming a median magnifi- 
cation of 10. The crossover luminosity is listed as "N/A" for the lensing survey because the clustering fluctuations are 
greater than the Poisson noise for all luminosities over which our model predicts sources should be detected. 



